Photon statistics and dynamics of Fluorescence Resonance Energy Transfer 



Andrew J. Berglund|] Andrew C. Doherty, and Hideo Mabuchi 
Norman Bridge Laboratory of Physics 12-33, California Institute of Technology, Pasadena, CA 91125 

We report high time-resolution measurements of photon statistics from pairs of dye molecules 
coupled by fluorescence resonance energy transfer (FRET). In addition to quantum- optical photon 
antibunching, we observe photon bunching on a timescale of several nanoseconds. We show by nu- 
merical simulation that configuration fluctuations in the coupled fluorophore system could account 
for minor deviations of our data from predictions of basic Forster theory. With further charac- 
terization we believe that FRET photon statistics could provide a unique tool for studying DNA 
mechanics on timescales from 10 -9 — 10 -3 s. 
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Fluorescence Resonance Energy Transfer (FRET) has 
become a widespread tool for probing molecular structure 
and dynamics. Recent demonstrations of single-molecule 
sensitivity in optical assays based on FRET jl], ^| have 
led to significant advances in our understanding of topics 
such as RNA folding and ribozyme function |3| . Interpre- 
tation of FRET data generally relies on a simple physical 
model involving near-field dipolc-dipolc interactions be- 
tween molecules, which was first proposed by Forster [Q. 
While some basic features of this model pertaining to 
steady-state solutions have previously been verified ex- 
perimentally, dynamical details have been largely inac- 
cessible. In this Letter, we report the use of a Hanbury- 
Brown Twiss apparatus to record photon statistics of the 
light emitted by FRET-coupled dye pairs with nanosec- 
ond resolution, and show that a careful comparison of our 
data with predictions of Forster theory supports the basic 
model but indicates a class of additional factors that must 
be considered. Our analysis suggests that conformational 
fluctuations of the substrate for the FRET-coupled dyes 
could be such a factor, which in turn points to the in- 
triguing possibility of utilizing FRET photon statistics 
for novel assays in DNA and protein mechanics. 

The FRET process involves non-radiative transfer of 
energy from a donor, which absorbs a photon of incident 
light, to an acceptor that is not directly coupled to the 
incident light. Detection of acceptor fluorescence is thus 
a simple indicator of FRET activity. A schematic energy- 
level diagram of molecular states is shown in Fig. [l| Un- 
der appropriate conditions of spectral overlap, Forster 
theory || predicts that the rate of energy transfer 
varies as Tp oc k 2 /R 6 , where n 2 depends on the ori- 
entation of the fluorescent species and R is the distance 
between them. For commonly-used organic dyes, the sen- 
sitivity of Tp to variations in R is greatest in the range 
of several nm; hence experimental determination of IV 
yields information on distance scales relevant to biolog- 
ical macromolecules. The R~ 6 distance dependence of 
Forster theory has been experimentally confirmed using 
'ruler' strands of DNA 0. 

Intensity correlation functions of the light emitted 
by FRET-coupled dyes pairs carries further information 
about the molecular physics of FRET. The second-order 



intensity correlation function (for a stationary process) 
gives the normalized time-average intensity Ik (t + t) of 
mode k at time t + r multiplied by the intensity in mode 
j at time t: 
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For j = k, this quantity is the autocorrelation of the 
intensity field j. Nonclassical photon statistics ||, for 

(2) (2) 

which 9jj (t) > g)j (0) for some r > 0, were first ob- 
served in an atomic beam f|. Since then, photon an- 
tibunching has been observed in a variety of systems 

(2) 

We have measured g)- k '(t) for individual FRET- 
coupled Cy3 and Cy5 dye molecules tethered to DNA 
ppf . We show that our data should be sensitive not 
only to the mean values of FRET parameters, but also 
to underlying molecular processes that perturb their val- 
ues on any timescale down to that of the radiative life- 
times. The techniques presented here thus provide po- 
tential experimental access to molecular dynamics that 
influence radiative level structures and couplings in the 
range 10-1000 ns, in a manner complementary to that 
of established techniques such as fluorescence correla- 
tion spectroscopy. Intriguing examples of such dynamics 
include photochemical processes, chemical shifts arising 
from changes in the local environment, single base-pair 
fluctuations in DNA secondary structure formation, and 
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FIG. 1: Energy-level diagram of molecular states relevant to 
the basic Forster model. Donor states on the left are coupled 
by FRET to acceptor states on the right. 
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FIG. 2: Schematic diagram of the apparatus. When making 
cross-correlation measurements, the 50/50 beam-splitters are 
removed to improve collection efficiency. Spectral filters and 
focusing optics at the APDs are not shown. 



conformational fluctuations on much shorter timescales 
than have previously been studied using FRET. Molec- 
ular dynamics simulations of nucleic acid mechanics are 
generally tractable only for integration times ~ 10 ns 
PI , so that experimental access to these timescales may 
provide fruitful contact between theory and experiment. 

A diagram of our apparatus appears in Fig. || It 
consists of confocal imaging optics coupled to Hanbury 
Brown-Twiss (HBT) detection channels Q. We fo- 
cus 140 /iW of 532 nm laser light between glass cover 
slips through a diffraction-limited microscope objective 
(Carl Zeiss). Fluorescence is collected by the same ob- 
jective and imaged onto a 100 ^m-diameter pinhole. A 
dichroic filter separates Stokes-shifted fluorescence light 
from scattered excitation light. A second dichroic filter 
separates donor fluorescence ( 570 nm) and acceptor flu- 
orescence ( 670 nm) into separate HBT channels each 
containing a 50/50 beam-splitter, spectral filters and 2 
avalanche photodiode single-photon counters (APDs). In 
each experiment, photon arrival times at one pair of de- 
tectors are recorded with sub-ns resolution by a time- 
interval analyzer (TIA). An electronic delay of 5 = 50 
ns is imposed in one channel to avoid small time-interval 
crosstalk in the TIA. 

We monitor fluorescence from dual-labelled DNA hair- 
pins in aqueous buffer at room temperature. The donor 
and acceptor are tethered at complementary positions, 
so that they exhibit a high FRET efficiency. In a typical 
experiment, we place 1/iL of 1 nM dye-labelled DNA so- 
lution between the cover slips. The axial position of the 
microscope objective is actively locked by a piezo-electric 
translator so that it is stable to < 100 nm for periods 
much longer than a typical experimental run (~ 6 hrs). 
We choose a low enough DNA concentration that there 
are no molecules in the imaging volume for a large frac- 
tion of the observation time. The count rate at the detec- 
tors then shows background light punctuated by bursts of 



fluorescence as individual DNA strands diffuse through 
the imaging volume. Since fluorescence from separate 
DNA strands is uncorrelated, the presence of multiple 
molecules in the imaging volume reduces the sharpness 
of features in the measured correlation functions. 

(2) 

To measure g y - k (t), we choose a pair of detectors 
and histogram the time interval between photon arrivals, 
keeping only data for which both detection channels were 
active (the TIA- limited channel dead time is 243 ns). We 
choose a threshold count rate for each channel (the max- 
imum expected count rate over a 1 ms interval, given 
the mean count rate over the ~ Is local measurement 
interval), and keep only those data for which at least one 
channel exceeds this threshold. In this way, we reduce 
the contribution from background light recorded when 
no molecule is present in the imaging volume. Mea- 
sured correlation functions, averaged over many molecu- 
lar transits, are shown in Fig. |3[ The precise time delay 
imposed by optical and electronic path differences (the 
t = point) is determined by measuring correlation func- 
tions of a pulsed LED. 

Interpretation of the experimental results proceeds 
from a straightforward model for Monte Carlo simula- 
tion of this and similar experiments. The donor and ac- 
ceptor are organic molecules attached to a complex sub- 
strate. Electronic excitation is followed by fast rotational 
and vibrational relaxation, so our model assumes negli- 
gible coherence between electronic states. We represent 
the donor and acceptor, labelled j — 1,2 respectively, as 
two-level emitters with basis states {|0)j, (see Fig. 
|l|) and lowering operator aj . Since we assume no coher- 
ent interactions, we write a master equation for the time 
evolution of the density operator of the system p with 
only incoherent (jump) terms: 
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dt 
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The 'jump operators' A m and associated (possibly time- 
dependent) rates Y m (t) represent the incoherent tran- 
sitions that can occur in the system. In our case we 
choose M = 5 possible transitions, but it is straight- 
forward to generalize the model to include more pro- 
cesses. These transitions are: direct excitation of the 
donor with rate r^; off-resonant excitation of the ac- 
ceptor with rate jT^, / < 1; donor (acceptor) spon- 
taneous emission with rate T\ (I^); and FRET energy 
transfer from donor to acceptor excited state with time- 
dependent rate Tp(t). A FRET transition is represented 
by Ts(t) = Tp(t), A5 = o"i <E> a\. The other jump op- 
erators are similarly defined. Eq. (|^) is equivalent to a 
linear system of rate equations for the ground and excited 
state populations of the donor and acceptor, and can thus 
be solved exactly when the time-variation of all param- 
eters is specified analytically. For stochastic parameter 
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FIG. 3: Measured correlation functions for high FRET- 
efficiency donor-acceptor pair fluorescence. Top: donor in- 
tensity autocorrelation (partially contaminated by inactive 
Cy5). Middle: acceptor intensity autocorrelation. Bottom: 
donor-acceptor cross-correlation, the average intensity in the 
acceptor channel, given a photon arrival in the donor channel 
at r = 0. See text for a discussion of the bunching at ~ ±5 
ns in the acceptor autocorrelation. 



FIG. 4: Monte Carlo calculation of correlation functions 
for a jump-process Forster transfer rate with correlation time 
tf = 7. Top: donor intensity autocorrelation. Middle: 
acceptor intensity autocorrelation. Bottom: donor-acceptor 
cross-correlation. The smooth curves are deterministic simu- 
lations with fixed Tp(t) = (Tf)- The parameters used were 



r x = r 2 = i, r" = 10, r 



0, / = 0.1. All rates (times) are 



in units of the laser excitation rate Tl (1/Tl). 



variation, we resort to numerical simulation. 

The measured autocorrelation functions in Fig. |3| show 
pronounced antibunching dips at r = 0. For a single 

(2) 

fmorophore, g)- ■ (0) = since two photons (in the same 
mode) can never be emitted simultaneously, i.e., in a 
time-interval r = 0. The observed value of (0) is a 
function only of the signal-to-noise ratio S at each detec- 
tor and the probability P(N) that N molecules are ob- 
served simultaneously. In order to understand the depth 
of the r = minimum in the autocorrelation functions, 
we make independent estimates of P(N) and S. Neglect- 
ing crosstalk between channels and assuming Poisson- 
distributed background, it can be shown from (Q) that 

N>1 Jy ' 

where Sj is the signal-to-noise ratio of HBT arm n, 
mode j. Assuming P(N) is a Poisson distribution with 
mean value (N), we determine the fraction of all pho- 
todetection events attributed to molecular fluorescence 
(i.e., exceeding the threshold count rate criterion de- 
scribed above). This fraction is J2n>i -^(-^0 from which 
we can solve for (N). For a typical run, we estimate 
(N) « 0.1 in this way. The probability that N > 2 
molecules are observed is therefore < 0.005 so that we 
are firmly in the single-molecule regime. We estimate the 



signal-to-noise ratios S)j by comparing the fluorescence 
count rate to the background count rate. A typical value 
is S ~ 5. From signal-to- noise and image volume occu- 

(2) 

pancy statistics, we expect g\ x (0) = 0.23 for the donor 

(2) 

and (0) = 0.32 for the acceptor autocorrelation. 

In our experiment, we see a large fraction (~ 60%) 
of low FRET-efficiency events, indicating a subpopula- 
tion of fluorophores exhibiting little or no FRET cou- 
pling. These events, which we attribute to acceptor pho- 
toblcaching, contribute a background that contaminates 

(2) 

the shape of g\ x (t). [In this experiment, we are lim- 
ited by the TIA to time-resolution on two detectors only 
and are thus unable to exclude FRET-inactive dye pairs 
when measuring g\ x (r).j Both the cross-correlation and 
acceptor autocorrelation depend on acceptor fluorescence 
events, and are therefore robust against a bare-donor sub- 
population. The depth of the observed r = feature in 
the acceptor autocorrelation is consistent with our in- 
dependent estimate based on signal-to-background and 
image volume occupancy statistics. 

In addition to photon antibunching on radiative 

(2) 

timescales, we see bunching (g\ 2 1 (t) > 1) at longer time 
intervals (~ 5 ns) in the acceptor autocorrelation of Fig. 
||. This bunching indicate clustering of acceptor events on 
the same timescale, most likely arising from fluctuations 
in TF(t). We expect that our fluorophores rotate with a 
characteristic time of ~ 250 ps, so rotational diffusion is 
an unlikely explanation for the observed correlations [ p"3"[ . 
Numerical simulations investigating the influence of rota- 
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tional diffusion on Tp(t) do not reproduce the features in 
our data. We rather suspect FRET 'intermittency', pos- 
sibly related to fast diffusion of tethered dye molecules 
and their propensity to stick to DNA (l4|, |l5|. Intersys- 
tem crossing and spectral diffusion for fluorophores such 
as Cy5 are known to exhibit longer timescales |lT| . 

We model FRET intermittency by allowing Tp(t) to 
jump between a high value Tp and a low value Tp (as 
perhaps when one or both dyes are stuck to the DNA) 
with a correlation time Tp. Numerical results are shown 
in Fig. ^, where we see antibunching followed by bunch- 
ing at Tp |2l|] . Deterministic simulations (smooth curves) 
with fixed (t) do not exhibit bunching in the accep- 
tor autocorrelation. Most calculated and observed fea- 
tures are consistent with intuition based on the four-level 
model. Under conditions of high FRET efficiency, donor 
emission rarely occurs. However, for a sufficiently strong 
driving field, the excitation rate is large compared to the 
acceptor emission rate, and donor absorption may occur 
when the acceptor is already in its excited state. In this 
"exciton blockade" situation, FRET cannot occur since 
the acceptor is already excited. Subsequent donor emis- 
sion is highly probable followed by acceptor emission a 
short time later. The conditional probability for accep- 
tor fluorescence is therefore enhanced by observation of 
a donor emission event, which is represented by a cusp 
at short positive times in the cross-correlation. The dip 
in the cross-correlation at negative r can be understood 
in a similar way. Observation of an acceptor photon de- 
terministically prepares the acceptor in its ground state. 
Since the acceptor is in its ground state, FRET occurs 
with high probability since any residual donor excitation 
is efficiently transferred to the acceptor. This depressed 
probability for donor fluorescence at short times after 
acceptor fluorescence is represented by a dip at negative 
r in the cross-correlation. Cross-correlations exhibiting 
other types of conditional statistics have been observed 
in cascaded multi-exciton emission from semiconductor 
quantum dots. |l9[] . 

In summary, we have measured FRET photon statis- 
tics and presented an intuitive model for interpretation 
of such experiments. Our data are in basic agreement 
with simple Forster theory, but the shape of the accep- 
tor autocorrelation is strongly suggestive of additional 
dynamics at ~ 5 ns timescales. Numerical simulations 
show that this inconsistency is resolved by the inclusion 
of stochastic variations in the donor-acceptor coupling 
strength. We have suggested a possible molecular mech- 
anism for these fluctuations, but further experiments are 
necessary to characterize the biochemical details. 

Our model suggests that thermally excited bending 
modes of dye-labelled, rod-like molecules should be vis- 



ible in g\ k (r) for appropriate motional amplitudes and 
timescales. We hope to exploit this dependence in study- 
ing the conformational dynamics of semi-rigid DNA and 
synthetic proteins. Furthermore, future experiments in- 
corporating direct excitation of the acceptor fluorophorc 
will allow absolute determination of model parameters. 
We believe these techniques can be developed into an 
important new optical single-molecule method for char- 
acterization of macromolecular dynamics on nanosecond 
(and longer) timescales. 
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